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Abstract 


Our interest is in the scaled joint distribution associated with ^-increasing subsequences for 
random involutions with a prescribed number of fixed points. We proceed by specifying in 
terms of correlation functions the same distribution for a Poissonized model in which both the 
number of symbols in the involution, and the number of fixed points, are random variables. From 
this, a de-Poissonization argument yields the scaled correlations and distribution function for 
the random involutions. These are found to coincide with the same quantities known in random 
matrix theory from the study of ensembles interpolating between the orthogonal and symplectic 
universality classes at the soft edge, the interpolation being due to a rank 1 perturbation. 
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1 Introduction 


To motivate our study of random involutions, we first recall the corresponding problem for 
random permutations, the solution of which is known. Let Sn denote the set of the IV! distinct 
permutations of {1,2,..., N}. For each ir e Sjy denote the image of the number i € {1,2,..., N} 
by 7T(*). A subsequence of image points 7r(*i), n(i 2 ), • • •, 7r(ij) where 1 < i\ < • • • < ij < N is 
said to be an increasing subsequence of length j if 7r(*i) < 7r(i 2 ) < ■ ■ ■ < 7 r(ij). More generally, 
we say there is a ^-increasing subsequence of length j if it contains k disjoint subsequences 

tt( 4°) < vr(4°) < • • • < 7r(4?), 1 < 4° < • • • < 4? ^ ^ = !> ■ • • > k ) 

with Yli=iji = 3- Being disjoint these subsequences contain no common member. For a given 
(k) 

ir, let Lr N (jt ) denote the maximum length of all the ^-increasing subsequences and define 

AjJV) =L Sv ) (tt) 4^V) :=0- ( L1 ) 


Note that 


> ■■■ > > 0. 


Consider an ensemble of permutations of {1,2,...,IV} in which each permutation is equally 
likely. The problem of interest is the computation of the scaled joint distribution of { (^)}j=i,...,i 

in the limit N —► oo. 

In the case l = 1 it was proved by Baik, Deift and Johansson [lj that 


lim Pr 

N— KX) 


A^ - 2 VN 

nv 6 


< s) = F 2 (s), 


( 1 . 2 ) 


where F 2 (s) is the scaled cumulative distribution of the largest eigenvalue for large random 
Hermitian matrices with complex elements (technically matrices from the Gaussian unitary 
ensemble (GUE)) [1011251. The latter is specified in terms of a Fredholm determinant according 
to 

F 2 (s) = det(l - K S0ft X(S}00) ), 

where xj is th e characteristic function of the interval J, while A' soft is the integral operator with 
kernel given in terms of Airy functions according to 


K soit (x,y) = 


Ai(x)Ai \y) — Ai(y)Ai'(x) 


f 


x-y 

i 

Ai(.x + t)Ai(y + t) dt. 


(1.3) 


(The superscript ‘soft’ indicates that the eigenvalue density pn)(s) is not strictly zero for any s, 
even though the origin is chosen in the neighbourhood of the largest eigenvalue and thus at the 
spectrum edge). The case of general l was solved in [?, 1191 !211 , where it was proved that 


lim Pr 

TV— XX) 


Aff - 2 VN 


N 1 / 6 


<s 


\ { $ - 2 VN 


AT/6 


< Sl^j = F 2 (s!,...,Sl). 


Here F 2 (s \,..., s;) is the scaled joint distribution of the l largest eigenvalues in the GUE. The 
latter is uniquely specified in terms of the /c-point correlation function 


p^(s 1 .»,)]«=! 


soft / 
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for the scaled eigenvalues at the soft edge of the GUE. 

It is the objective of this study to calculate the analogous joint probability for involutions 
with a prescribed number of fixed points. One recalls that an involution is a permutation 7 r with 
the additional property that 7r 2 = /, where I denotes the identity permutation. Involutions 
must consist entirely of two cycles and fixed points. Thus for an involution of {1,2,..., N}, if 
there are n two cycles, there must be m = N — 2n fixed points. As emphasized in f2j, a random 
involution with a prescribed number of fixed points can be generated geometrically by marking 
n points (n < [N/ 2]) in the unit square below the diagonal y = x uniformly at random, marking 
the images of these points under reflection about y = x, and marking N — 2n points uniformly 
at random on the diagonal. Equivalently the unit square can first be divided into a N x N 
integer grid, and the points marked at random on the lattice sites below and on the diagonal 
according to the above prescription, with the additional constraint that no two points are in 
the same row or column. Either way, projecting the points onto the x-axis gives a sequence of 
x co-ordinates x± < X 2 < ■ ■ ■ < x/v while projecting them onto the y -axis gives a sequence of 
y co-ordinates y± < y 2 < ■ ■ ■ < yN■ Each point will then have a co-ordinate (xi,y n ^) with the 
property that 7 r(z) = i for the point on the diagonal, and that (x^),?/*) is the point reflected 
in the diagonal otherwise. Hence ir defines an involution with a prescribed number of fixed 
points, and furthermore the involutions are generated at random with uniform probability by 
this procedure. 

The quantity L^\n) admits an interpretation in the above setting. Connect points by 
segments which always have positive slope to form a continuous path, which is said to be 
right/diagonal (rd). Define the length of this path as the number of points it contains, and 
denote it by #rd. Similarly, let (RD) fc denote the set of all k disjoint rd lattice paths, formed 
from amongst the points with a path being disjoint if it contains no common points, and for 
(rd) fc £ (RD) fc let #(rd) fc denote the number of lattice points. Then by considering recurrences 
satisfied by the various quantities one can show .22 a 

L^tt) = max ^ #(rd) fc (1.4) 

(rd) fc e(RD) fc 

(for the recurrences satisfied by the left hand side, see |T2] ( Appendix A]). The equation (11.41) 
holds for points corresponding to a general permutation 7 r. In the case that 7r is an involution, it 
is easy to see that the set of lattice paths (RD) fc can be restricted to those which contain points 
on or below the diagonal (see Figure Q] for an example). 

The explicit form of the scaled distribution of =: A n)m for random involutions with n 
2-cycles and m fixed points, and thus the solution to our problem in the case l = 1, is already 
known 3]- Thus introduce the scaling variable w by the requirement that 

m = [y/2n - 2w{2 n) 1/3 ], (1.5) 

where [•] denotes the integer part. It is proved in [3| that with w fixed 

jv 1 !So Pr (^ r W^ £s ) =F "( s; ’" )l (L6) 

where the distribution F a (s\w) is specified in terms of a certain Riemann-Hilbert problem 
related to the Painleve II equation with special monodromy data. This distribution has the 
property that 

F a (s-0) = F 1 (s), lim F^(s-w) = F 4 (s), 

W —>00 
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Figure 1: Nine points in the square symmetric about the diagonal which correspond to the 
involution (12)(3)(46)(59)(7)(8). An rd path realizing has been drawn. 


where Fi(s), F^s) are the scaled cumulative distributions for the largest eigenvalue in the GOE 
of random real symmetic matrices, and the GSE of Hermitian matrices with real quaternion 
elements respectively. 

Here we seek the joint distribution of {An'L(7r)}j=i,...,z, scaled as in (11,51) . (11.(11) for general l. 
Regarding the joint distribution of all the A j = 1,..., N as specifying a point process, 
it is generally true that the distribution of the l right-most points is fully determined by the 
corresponding correlation functions. The specification of the correlations for a point process 
implying the distribution F a (s',w) for the right-most point has recently been given |11| . This 
was found in the study of a closely related geometrical model to that corresponding to random 
involutions and their increasing subsequences. Thus consider an M x M integer grid. Associate 
with each lattice site on or below the diagonal a continuous exponential variable 

p r (xij e[y,y + dy}) = e~ y dy, i < j, 

Pr(x M € [y,y + dy]) = ^ 2 ^ e ~ (1 ~ A)y/2 dy, (1.7) 

and impose the symmetry constraint that x l .j = xy,; for i > j. The quantities m, with 

the lattice points counted according to their weightings xyj, are well defined. Setting then 
(k—1) (— 1 ) 

Xk = L m — L m , L m := 0 we have x\ > X 2 > ■ ■ ■ > xm > 0 and furthermore the joint 
distribution of these variables is proportional to fl Section 3] [T2J Prop. 4] 

e -Ef =1 ^/2 e A E r =1 (-ip-G,/2 Yl (Xi-xj). (1.8) 

This same p.d.f. occurs in random matrix theory. Thus let X be a 2n x 2n antisymmetric 
complex Gaussian matrix (independent entries distributed as N[0,1] + iN[0,1]), and let x be a 
2n x 1 complex Gaussian vector with entries distributed as N[0, l/\/2] +?'N[0, l/y/2]. Then we 
know from Thm. 3] that Y = X^X + bxx 1 has eigenvalue p.d.f. (11.81) with A = 1 — 2/6, 
M = 2 n. 
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For the p.d.f. m, the scaled correlation functions with 

Xj = 4 M + 2(2M) 1 / 3 X j , A = u/2{2M ) l/3 , 


(1.9) 


and M —> oo have been computed T1 (actually the symbol a is used instead of u in m, but 
we use a for another purpose below). On the other hand, a result of Baik [Tl gives 


lim Pr d«h-w<2M)-/» 4M < „) = F * (s . w y 

M—>oo V 2 2M l 1 / 3 / y ; 


2(2M) 1 / 3 

The scalings m is precisely that in (ITTol) with the identification 


u 


w = — - 


( 1 . 10 ) 


( 1 . 11 ) 


so this distribution is fully determined by the correlations computed in m • We therefore expect 
that for random involutions, the generalization of c@ is 


lim Pr 

N—>oo 


X ( n) n -2 VN „ XSm-2 VN 


N 1 / 6 


< Si,..., 


N 1 / 6 


< Si= F 0 (si,... ,si;w ), 


( 1 . 12 ) 


where F 0 (si,..., sf, w) is the joint distribution of the scaled l right-most points in the process 
specified by tea, scaled as in m and with u and w related by cm As revised below, it is 
generally true that the latter distribution is fully determined by the scaled correlation functions 
P S (kf ed - According to jTT] these have the explicit form 


pf k f ed (Xi, • • •, X k -u) = qdet .V,)],, ,. 


(1.13) 


where with sgn(u) = 1, —1, 0 for u>0, u<0, u = 0 respectively, / is a 2 x 2 matrix with 
entries 


f U (X,Y) 
f 22 {X, Y) 


f 22 (Y,X) 

1 
2 

,-x 


K soit {X,Y)~ -A [ X e u{x ~ t)/2 K soit {t,Y)dt 
2 oY .1—00 


u 

’4 


/ °o a 

ds X-K soh (s,t), 


1 (u d 


= —e 


4 V 2 dX 

u\X— 


U 


d 


x 


= 4li + avlli + aF)i / 

f\X,Y) 


K 


soft 


lY 




cX 


y|/2 sgn(X - Y) - { [ e^-P/'^K^iX^dt 
^ J — OO 

e u(X-t)/2 K soft(Y_ t) dt^. 


(1.14) 


and the qdet operation is defined in terms of the more familiar Pfaffian by the general formula 

^2 n = 1 n 


qdet A = Pf (AZ. 1 


2 n h 


0 -1 

1 0 


valid for all 2 n x 2 n matrices A with the self dual property A = Z 2 y A T Z 2 n - Our main task then 
is to show that the same scaled correlations determine the joint distribution of {An'L(7r)} J=1) 

We will see that doing this indeed allows cm to be validated, giving us the following limit 
theorem, which is our main result. 
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Theorem 1. Consider an involution on N = 2 n + m symbols, consisting of m fixed points and 

(k) 

n 2-cycles, chosen at random with uniform distribution on the set of such involutions. Let A „.J n 

(k) 

be specified in terms of the maximum length of the k-increasing subsequences L K N according to 
<Z3>. With m related to n by the limit formula EH holds. 

2 Random Generalized Involutions 

2.1 Strategy 

We are guided in our approach by previous studies on the calculation of distribution functions 
relating to increasing subsequences for random permutations mm, and also previous studies 
on the calculation of the distribution of the maximum increasing subsequence for random invo¬ 
lutions HE]. For both random permutations and random involutions, all the studies proceed 
by Poissonizing the ensemble. In the case of random permutations, the number of symbols N is 
itself taken as a variable which occurs with probability e~ z z N /N\. For random involutions the 
number of two cycles n and number of fixed points m are separately taken as random variables, 
occurring with probability e~ Zl ~ Z2 z'fz^ l /n\m\. One then seeks to calculate the distribution func¬ 
tion of interest in the Poissonized ensemble, and to compute the z —► oo (z\ , zi —■ ► oo) scaled 
limit. A de-Poissonization argument m gives that the latter is equivalent to the N —> oo 
(to, m —* oo) limit in the original ensemble. 

As the distribution functions are controlled by the correlation functions P(k), one would like 
to first calculate the Poissonized correlation functions. However we know from studies of random 
permutations that a direct calculation is not practical. Instead, as was shown by Johansson m, 
progress can be made by constructing the Poissonization as a limiting case of a generalization 
of the geometrical viewpoint of the increasing subsequence problem. In one such generalization, 
giving rise to the so called Meixner ensemble m, each lattice point of an M x M grid carries 
a non-negative integer variable chosen from the geometric distribution with parameter q. For 
this model pn c ) can be computed in terms of certain orthogonal polynomials. Taking the limit 
M —> oo with q = Q/M 2 gives p^ for the Poissonized version of the original model. Thus we 
must first calculate p^ for the version of the Meixner ensemble relevant to involutions. 

2.2 The joint p.d.f., correlations and distribution functions 

According the above strategy, our first task is to introduce a generalization of the geometrical 
model corresponding to random involutions and their increasing subsequences. For this we 
consider an M x M grid, and associate with each lattice site on or below the diagonal a non¬ 
negative integer variable 


Pr (xij = k) = (1 - q)q k , i < j, 

Pr (x iti = k) = (1 - ^aq)(aq) k/2 . (2.1) 

For i > j we impose the symmetry constraint = Xj j. Notice that the symmetry constraints 
of this model are the same as that for the geometrical model of random involutions specified in 
the Introduction, and notice too the similarity with the model defined by 0) and surrounding 
text. As with the latter model, the quantities Q, which can be regarded as a sequence of last 
passage times, are well defined. With A& := L := 0, we know from [TJ Section 
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3 ], m Prop. 1] that the joint distribution of the latter variables have the explicit form 


P M (A;?) = (1 - y/aq ) M (1 - ^(M-p/2^ J] 


1 <j<l<M 


\j — Xl + l — j 

—J ' 


In terms of hj := Xj + M — j this reads 


where 


and oo > h\ > > ■ ■ ■ > h,M > 0. Note the similarity between El and m- 

Let us introduce the symmetrized joint p.d.f. by 

^sym(^'l) • • ■ ) Hm) = P(hn{l) j ■ . . 1 h/u(M))Xh /J ( 1 )>--->h M ( N ) j 
A 



(2.2) 

(hj - hi), 

1 

(2.3) 

M -1 

n (wo 

i=i 

(2.4) 


(2.5) 


where \T = 1 for T true, and \T = 0 otherwise. The Appoint correlation pn.\ is then given by 


P(k ) {hi j • • • j hfX) 


1 

(M - k)\ 


oo 

Y P(hi, ■ ■ ■, Hm)- 

hk+i,—,hM=0 


( 2 . 6 ) 


The correlations El can be used to compute the distribution functions Pr(/ii < a\,... ,hi < 
ai ) for the joint p.d.f. (12.51) . To see this, let ao = oo, let a\ > <22 > • • • > a/ be non-zero integers, 
and put Ij = (aj,aj-i) where (aj, dj-\) denotes all integers between (but not including) aj and 
dj- 1 - Let EM({(n r , I r )}r=i,...,l) denote the probability that n r of the coordinates {hi}i=\,...,M 
are in I r (r = 1,..., l). Then as a consequence of the definitions we know from m eq. (3.41)] 
that with 

r 

L z := {(m,... ,n z ) G Z> 0 : Y rij < r - 1 (r = 1,..., l)}, 

5 = 1 

one has 

Pr(/ii < ai,... ,hi < a{) = Y E M ({(n r , I r )}r=i,...,l)- (2-7) 

(ni,...,n()eLi 

Furthermore, again from the definitions, it is easy to see that 


-^Af ({(^r> 7r)}r= 1 ,..., Z ) 


( — \)J2 r =l nr Q^j=l n i 
n\ \ ■ ■ ■ n z ! 9£ ni • • • <9£ n * 


M k 


(nO-E«fof)) 


Psym 


5i=-"=5fc =1 


( 2 . 8 ) 


(?) • (?) 

where Xj T = 1 if hj G I r , Xj r = 0 otherwise, and that the average herein is given in terms of 
the correlations El according to 


M k 

(n(i-M>)) 


Psym 


M , 1 v p 00 p k 

1 + E^r U(Y^X? r )p( P )(hu...,h p ). (2.9) 

p= 1 hi,...,h p =0 i=l r=l 
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2.3 Quaternion determinant expression for the /c-point correlation 

We know from DU Eq. (3.1)] that with P given by (12.311 . and M even (12.51) can be written 

M 

-fsym(^i, • • • ,^m) = a)q^i =1 h ^ 2 JJ(/r,- — hi)Pi[e(hj, /iz)]j,j=i,2,-,M> (2-10) 

j>i 

where 

e(x, y ) = a^ -x l/ 2 sgn (y — x ) (2-11) 

(a similar formula can be given in the case that M odd, but for definiteness we will proceed with 
the assumption that M is even). This structure is familiar in random matrix theory, and in the 
theory of random measures on partitions relating to increasing subsequences m- It is known 
mm that in general the corresponding /c-point correlations have the quaternion determinant 
form 

P(k)(hi, ■ ■ • , h k ) = qdet[f(hj, /i/)]j,z=i,2,- (2.12) 

where f(x,y ) is the 2x2 matrix 


f(x,y) 


S(x,y ) I(x,y) 
D(x,y ) S(y,x) 


(2.13) 


According to the general formalism, to specify the matrix elements in (ESP we must introduce 
skew orthogonal polynomials R n (x ) satisfying 


(^2m(^)i 7?2ra+l(?/)) — {R<2n+l(x')i R2m(y)) — Tm&m,ni 


(^2m(^)i R^niy)') — 0) (^2m+l (®) j ^2n+l (?/)) — 0, 


(2.14) 


where 


( f(x),g(y )) = -^2^q {x+v)/2 e(x,y){f(x)g{y) - f(y)g(x)) 


x=0 y =0 


J2<i v/2 a~ y/2 f(y) E ^ /2 « x/ 

y=0 


One then has 


where 


(M/2)—1 


,2 9(x) ~^2q y/2 a y/2 g(y) E Q x/2 a x/2 f(x). 

x—y -\-1 y=0 ai=7/+l 

(2.15) 

(2.16) 

(2.17) 

(2.18) 


f(x,y ) = 

(M/2)—1 


5'(^,y) 7(x,y) 
^(x,y) S(y,a:) 


S'(x,y) = g y/2 — [$^(x)i2 2 j+i(y) - $2j+i(®)-R2j(j/)], 


r 9 

j=o 3 


I (x,y) = ~ E — [ $ 2i(x)4> 2 /+i(2/) - 4> 2j+ i(x)4>2j(y)]+e(x,y) 


j=0 


and 


(M/2)—1 


L»(x,y) = g^/ 2 £ -[i2 2j -(x)i2 2i+ i(2/)-i22j + i(s)i?2i(y)] (2.19) 


— r 7 
j=o 3 






with 


( 2 . 20 ) 


®j( x ) = ^2 e (y> x )(i y/2R j(y)- 

y =o 

2.4 Skew orthogonal polynomials 

The explicit forms of polynomials with the skew orthogonality property with respect to 

the skew product (12.151) is required. These can be calculated from general determinant formulas 
for skew orthogonal polynomials. 

Lemma 1. Let ( , ) be a general skew symmetric product. Let {Rj(x)}j=oA,... be the correspond¬ 
ing skew orthogonal polynomials which thus satisfy [2.Iff) - Let {Cj(x)}j = 0 , 1 ,... be any family of 
monic polynomials. Assuming V n and £ n as specified by \2.24\) and 1051) below are non-zero 
we have 


R2n(x) = P; 



C 2 n (x) 

j2n 2n—l 

j2n 2n—2 

j2n 0 

Vn 1 

C 2n -i(x) 

j2n—l 2n—l 

j2n—l 2n—2 

j2n—l 0 


C 0 (x) 

jO 2n—l 

jO 2n—2 

• j 00 


( 2 . 21 ) 


and 


R2n+1 (x) — E, 


-1 


j2n+l 2n+l 2n —1 j2n-\-l 2n—2 

j2n 2n+l J^ n J^ n ^ n ~2 


j2n+l 0 
j2n 0 


where 


jo 2n+i c 0 (x) J° 2n_1 J° 2n ” 2 

+c n R 2n (x), 

J mn = (C m (x),C n (y)), 


J { 


o o 


p». = 


f — - 


j2n—l 2n—l 

j2n—l 2n—2 

j2n —1 0 

j2n—2 2n —1 

j2n—2 2n—2 

j2n—2 0 

jO 2n—l 

jO 2n—2 

j 00 

2n+l 

j2n 2n—\ 

j2n 0 

j2n—l 2n+l 

j2n—l 2n—l 

j2n—l 0 

jO 2n+l 

jO 2n—l 

.. J0 0 


( 2 . 22 ) 

(2.23) 

(2.24) 


(2.25) 


and the c n are arbitrary constants. 

Proof. We note the skew orthogonality properties are equivalent to the requirements 

(Cj(x),R 2 n(y)} = (Cj(x),R 2 n+i(y)) = o, j <2n — 1 (2.26) 

and 


( C 2n (x),R 2n+ i(y )) / 0. 


(2.27) 
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The property is immediate from the structure of (non and non . while (1071) follows 

from the assumption that T> n+ \ is non-zero. Furthermore, by inspection both (toil) and (1071) 
give monic polynomials. □ 

The structure of (1071) suggests choosing {Cj(x)}j=o j i ! „, as particular Meixner polynomials. 
We recall EU the Meixner polynomials are defined as 


M n (x;c,q)= 2 F\ 


—n, —x 


and satisfy 


^ ~T qXMm ^ x 'q)M n (x-c,q) 


1 - 


q n n\ 


x=0 


(c)n(i - q) c 


with (c) n = T(c + n)/T(c). We define monic polynomials C n (x ) as 


n\q n 


C n {x) = j- — n M n (x\ 1, q) = x n + • 
From (12.291) the orthogonality relation for the C n {x) is 


x=0 


where 


^ ) q Cm{F)Cn{x) — hnSm,nj 
( n\) 2 q n 


h n. — 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


(1 - q ) 2n + 1 ‘ 

Using the symmetry M n {x\^,q) = M x (ri] 7 , q), (12.291) also implies the completeness relation 


(2.32) 


^ ^ , Cn(%)C n (y) — q Sx,y> 2-j V ^ Z> 0 . 


n=0 


h 


With {Cj(x)} so specified we seek to compute (12.231) . 
Proposition 2. For m < n, 

jmn = jnm = ( Cm (x),C n (y)) = a m b n , 

where 


bn = n 


y/q( 1 - y/aq) 

a- yfq V(! - o)(V<x - Vo) 
, 1 f Vq(V®-yq) \ 

'i-vm v(i-?)(i- v°®) 


Proof. Set 


(2.33) 


(2.34) 


(2.35) 

(2.36) 


F(y;a) = q y/2 a y/2 ^ q x/2 a x/2 C n (x). 

x=y+l 

It follows from the fact that C n (x) is a polynomial of degree n in x that F(y ; a) is a polynomial 
of degree n in y (thus y can therefore be regarded as a continuous variable). We can therefore 
write 

n 

F{y\a) = ^2 K n jCj(y). (2.37) 

3 =0 
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To calculate the n n j we first multiply both sides by q y / 2 a y / 2 and subtract the same equation 
with y replaced by y — 1 to obtain 

C n (y) = ~F{y\a) + q~ 1/2 a~ 1/2 F(y - 1; a). (2.38) 


To proceed further we note, as can be verified from the definitions (IOS1) and (EH, that 


Q 77,2 (j2 

xC n (x - 1) = C n+ i(x) -— (2n + 1 )C n {x) + - —rC n -i(x), 


q~ 1 


(q-iy 


n < \ n i \ nc t + n + q n , . n 2 q 
xC n (x ) = C n +i(x) --- C n (x) + -- — 


9-1 


(9 - i) s 


(2.39) 


(the second of these is the three term recurrence, the general structure of which holds for all 
sequences of orthogonal polynomials). Thus with (12.371) substituted in (12.381) . we can multiply 
both sides of the equation by y and use the equations to get an equation involving only the 
linearly independent functions C n+ i(y), C n (y),, Co(y) in the variable y. Equating coefficients 
of these functions shows 


\fq_ _ _ n\ Jag ( y/q(l ~ y/^Q) \ 

y/a-y/q' Knk k\ (y/a - y/q) 2 \(1 - q)(y/a - y/q)) 


(2.40) 


We will now make use of jOH> and (IOU1) to evaluate dosii . Substituting (1071) in (IOH1) . 
multiplying both sides by q v C rn (y ) and summing over y, making use of (12.311) on the right hand 
side, we see that 


q y/2 a y/2 C m (y) ^ q x/2 a x/2 C n (x) = 

y =0 x=y -\-1 

Recalling and making use of the explicit formulas (1071) and sim . (lorii follows. □ 

In general CD Prop. 6], if the skew product has a factorization (12.341) for rnonic polynomials 
C n (x), n = 0,1,..., the corresponding monic skew orthogonal polynomials Rk(x) can be written 
as a series in {C n (x)} n= o,i,...,fc for explicit coefficients involving the afis and bfi s, while r n = 
« 2 n^ 2 n+i- These facts can be seen from the determinant formulas (12.211) and (12.221) . We thus 
have the following result. 

Proposition 3. Let C n (x ) be specified by \2.ti(A) . The monic skew orthogonal polynomials with 
respect to nnsu are given in terms of these polynomials by 



hrn Knrn ; 1X1 XI, 
0, m > n. 


(2.41) 


R2n(x ) 


V®-yq (2n)\ q n-k-(l/2) 

1 — y/aq (2k + 1)! (1 — g) 2n - 2fc_1 2k+1 


R 2n+ i(x) = C 2n +i(x ) - -y- n + l)C 2n (x) 

1 - y/aq 1 -q 

while the corresponding normalization has the explicit value 

(2n)\(2n + l)\y/Fcq 2n+ W 


r n = 


(l-q)^+\l-y/aqf 


(2.42) 

(2.43) 


(2.44) 
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Proof. These results follow from the general formulas of EH. The only point which requires 
clarification is that we have chosen for the arbitrary constant c n in 


c „ = A4n^vT (2n+1) , 

1 - y/aq 1 - q 

which gives the simplest expression for R. 2 n+i(x) in the basis {C'j(.r)}. □ 

We note that following the procedure used in m , the triangular structures (TP 21 ) and ll!P31l 
can be inverted to give 


C 2 n(x) = R. 2 n(x) + 


7 2n (2 n)\q n 


(1 - q) 2n 

n ~ l (1 -q) 2k 
7 '* ) jf —■* 'y 2k (2k)lq k 


E 


n— 1 


R 2 k(x) + ^2 


(1 -q) 


2k+l 


—' 7 2fc + 1 (2fc + 1 )!g fc +( 1 / 2 ) 


R2k+l(x) 


(2.45) 


and 


gl/2 ^ ^ ^ ^ ^ ^,2n+l (2n -j- l^]g n +(l/ 2 ) 


C 2 ri+ 1 (x) = 7 _—( 2 n + l)ii 2 ri .(x)+ (1 _ g) 2n+i 


1 \ ? \ (1 — q) 2k n ' ^ _ q^2k+i 

1 ~ 7 1 J ^ + ^ 7 2fc+1 (2fc + 1)!^+(V2) ^+1^ 

—0 k —0 


where 


7 = 


« - y/q 

1 - yfaq ' 


(2.46) 

(2.47) 


2.5 A summation formula 


Let us write 

n 

C n (x) =y1/3 nj Rj(x), (3 nn = 1. (2.48) 

3=0 

We see from the explicit formulas m . m that the coefficients have the factorization 
property 

0nj = Otnlj (2.49) 

for certain a n , 7 j. According to |23l Eq. (3.42)], under such a circumstance the matrix element 
S(x,y ) can be summed. 


Proposition 4. 

S(x,y) = 


The expression can be simplified to read 

(x+y) /2 1 Cm(x)Cm-i{U) - CM-l(x)CM{y) y /2 1 1 ~ Q 

hM- 1 a;-y h^-2)/2 M 7? 1/2 

X ^m- 2 (x) - g" 72 ^ 2 ^-!^)) (CW(y) - i?M(y)) • (2.50) 
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Proof. We know from general formulas m that with {(3 n j}j=o,...,n the lower triangular 
transition matrix for the change of variables from the monic skew orthogonal polynomials 
{i?j(a;)}j = o,i,... to the monic orthogonal polynomials {Cj(x)}j= 0 , 1 .... as implied by (12.481) . the 
quantity (ingnii permits the expansions 


$2fc-i(a0 = -q x/2 J 2 


is=2k—2 
oo 


$2k-2(x) = q x/2 X 


C v {x) 

hi? 

C u {x) 

h,. 


Pv 2k—2 Tk- 1 , 


Pv 2k—1 Tk- 1 - 


(2.51) 


u=2k—l 

Further, we know that these formulas substituted into dz> give 

S(x, V) = ,( ^ f) £ (2.52) 

!2=0 ^ I2=M fc=0 " 

According to the Christoffel-Darboux formula from the theory of orthogonal polynomials, the 
first summation gives the first term on the right hand side of Moreover, making use of 

the factorization (127731) shows 


u-M (v-M)/2 , 

X PukRkiy) = ^_ q y-M - Rm(v)) ■ 

After further use of (12.511) the second term of (127501) results. □ 

3 Poissonized Random Involutions 

3.1 Poissonization 

Essential to the rationale underlying the introduction of the geometrical model in 12.21 is its 
reduction in an appropriate limit to the Poissonization of random involutions. To see how this 
comes about, the procedure for the Poissonization of random involutions given in [2J must be 
revised. 

Consider involutions of {1, 2,... , JV}, and specify the number of fixed points therein to be 
m. Catalogue the involutions according to (A^ =: A k)k=i,...,N where is specified by (11.11) . 
The Robinson-Schensted-Knuth correspondence (see e.g. El) tells us that there is a bijection 
between all involutions with a given value of A = (Ai, A 2 , • • •, Ajv) and standard tableaux of 
shape A and content N. Furthermore, the number of fixed points m is equal to l) -7 " _1 A_j, 

and so restricts the permissible A. 

Let the number of standard tableaux of shape A and content N be denoted f\. With i(X) 
denoting the number of non-zero parts of A, and 

U(A):= n (A, — A, + (- j), W r ( A) = n . (3.1) 

1 <j<l<p j=li¬ 

lt is known that El 

fx = N\V l{x) (X)W l{x) . (3.2) 
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Further, let t njm denote the number of involutions of {1,2,... ,N} with m fixed points and n 
2-cycles. We have 


t 


n,m — 




Aen 


(oo) 

2n-\-m,m 


(2 n + m)\ 
2 n n\m\ 


(3.3) 


where 


(oo) 

2 n+m,m 


is the set of partitions with 


i{\) 

A j := | A| = 2 n + m, 

3 =1 


l(X) 

E(-irA 


j=l 


= m. 


(3.4) 


In terms of these quantities the probability that such an involution corresponds to a standard 
tableau of shape A is then given by 

Qn,m(^) = T X Ap n(°°) 

Lji jji L 2 n+m,m 

And the Poissonized form of this, in which n and m are regarded as random variables from 
distinct Poisson distributions, is 


Q(*i.*a)(A) 


zi—Z2 z 'n z ™' 1 


c>-z\-Z2 y n m 

e Zi %2 


fx 

n,m 

2 n 

( 2 n + m)! 


/a, 


(3.5) 


where the second equality follows upon using (EHl) . and n, m are related to A by JO). 

For future applications it is convenient to write this in terms of ( N,m ) rather than (n, m). 
Setting 

2 i = Q/2, z 2 = \[Qa (3.6) 

and with |A| = N defining 


p(Q,a)(A) 


= g -V^Q-Q/2 


a m/2Q N /2 


N\ 


fx, 


(3.7) 


we see that 

q(Q/ 2 ,(c*< 2) 1/2 )^ _ p(Q,aOm. 


(3.8) 


Also for future application, we note from the fact that P^ ,c 0 is a probability distribution that 
one has the summation 



N =0 


— gV aQ+Q/2 


(3.9) 


where 

t^:= (3-10) 

A:|A|=iV 


(the case a = 1 is given in (201 §5.1.4, eq. (42)]). 

Consider now a function ( 7 (A) of the parts of a partition A. A concrete example of future use 
is 

OO 

9W = n( 1 + u ( x j-j))> ( 3 - n ) 

3 = 1 
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where u(x) is any function which vanishes for x < 0. For such functions we define the Poissonized 
average as 

(5 (A)) (Q ’ q) =^ 5 (A)P {Q ’ q) (A). (3.12) 

A 

We want to show that this average results as a limit of the model of m 

Proposition 5. Let g( A) be a function of the parts of a partition X which furthermore satisfies 
the bound 

m\ (3-i3) 

for some c > 0. With Pm(A;(/) given by \2.X\) we have 

lim E g(X)P M (X;Q/M 2 ) = (g(X))^ a \ (3.14) 

M—>oo *—* 

Proof. We essentially follow Johansson 1X9. who proved the analogous result for Poissonized 
permutations. Using the explicit formula Q , we see 


E g{X)P M {X-q) 

X1(X)<M 


M M 


(1 -^q) M (l-q) M{M - 1)/2 E 




n 


1 <j<k<M 


X:i(X)<M j =1 j =1 

Xj ~ Xi + l — j 

—j 


M—l 


(1-V^) M (1-9) M(M_1)/2 E q N/2 * m/2 g(X)v M (X) n \- 

X-.l(X)<M j =1 J ' 


(3.15) 


With denoting the set of partitions A with the properties (13.41) but constrained so that 

K(A) < M, recalling (13.11) we see 


M—1 


E g( a)^m(a) [] - 


Aeo 


(M) 


3 = 1 


Aeo 


^ 9(A) w Mi a) n * 


e »w4n (Aj+M ' J)! 


Aeo 


(°°) 

N,m 


3 = 1 


(M-j)! 


(3.16) 


From the assumed bound and the further bound 


(A + 

M x i(M-j)\ ~ 


(a consequence of Stirling’s formula), we see from (13.161) that with q replaced by Q/M 2 in (13.151) . 
the sum over A : £(X) < M for M —> oo in the latter is itself bounded by 


^ Q N / 2 a m / 2 ^ 

jvT tN,m ’ 


Q = Qc 2 . 


15 







But according to (THUll and El this is summed as the r.h.s. of El with Q replaced by Q. 

Thus we have a uniform bound in N for the sum in EH with q replaced by Q/M 2 . The limit 

M —> oo can therefore be taken term-by-term. Doing so gives EH- □ 

We remark that eh can be written 

(9(A))<«’“> = (3.17) 

N =0 

<aW)P ■■= -JSi E s(\)a m,2 h- (3.18) 

A: | A | =iV 

Some insight into the role of a can be obtained by computing the mean number of fixed points 
in the ensemble EH - From EH, the identity El , the constraint involving m in El, and 
the explicit form El it is immediate that 


3 =1 


This in eh and use of El implies 


tW 

E<-F-A 

3 =1 


(“> r tN -1 N 
= Ja , s N. 


N 


t 


(«) 

N 


Recalling the explicit formula El , the sum EH can be estimated to give tffi/tff_ l ~ VN 
as IV —*• oo (a > 0). Hence 


(“) 

rsj 

N N —>oo 


V aN, 


E(-i rA 

3 = 1 

which tells us that in the ensemble eh the mean number of fixed points is necessarily pro¬ 
portional to y/N, with proportionality constant y/a. 


3.2 A Poissonized average and limiting correlations 

According to the equations El-El, to calculate the distribution function Pr(/ii < ai,..., hi < 
ai ), it is sufficient to compute the average 

M k 

1=1 r= 1 

Let us suppose each end point aj in the {I r } is a positive integer when measured from M. We 
indicate this by the replacement I r i—> M + I r . The quantity being averaged in eh is now of 
the form eh , and eh tells us that 


M 


(0 


^(nb-E&v 


1=1 


r =1 


Psym 


q=Q/M 2 


nd-E 

1=1 r=1 


0 a,Q ) 


(3.20) 
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We know from the proof of Proposition 0 that the convergence is uniform in M , and can 
therefore be taken term-by-term in (EH to give 


II f 1 ~~ ^Z&XXi-ieir 


1=1 r=\ 

oo 


(a,Q) 


/ , oo P k 

l + Y,~^r n ( Y ^rXj!) P {p) (h x + M, ...,h p + M) 


p =i 


pi 


hi,...,h p =0 i=l r =1 


q=Q/M 2 


(3.21) 


Moreover, the uniform convergence tells us that the limit M —* oo can be taken term-by-term in 
the summations of the quantities specifying p^ (recall (12.1711 (12.201) 1. To calculate these limits 
requires an appropriate asymptotic formula relating to C n (x). 


Proposition 6. For M —»• oo we have 
C n (x + 


n\q( n M -r)/2 


J—n+M+x (2-^/Q) 


q=Q/M 2 (1 - q)' 

uniformly in n G Z>o- 

Proof. Following CHI, we use the generating function of the Meixner polynomials m 


(3.22) 


1--) (l-0 _ 

9/ 


to derive the integral representation 


oo / X 

Y —r M n{x\c,q)t r ‘ 


(3.23) 


M n (x-c,q ) = 


n! 


2iri(c) r 


dt t 


n =0 


—n—1 


i-- (i - 0“ 


(3.24) 


where the path of integration encloses the origin anticlockwise. It follows from dSHH) that 
C n {x + M)\ q= Q / M 2 

dt t - n “ 1 1 1 M M ^-*--^-1 


n\q n 


2iri(q — 1) 

n \q(n—M—x)/2 ^_-y^~n+M+x 


1 - - 
q 


(i - q) r 


2vr 


(1 ~t) 
dfl Z ~ n + M + X 


Mz 


x+M 


1 - 


VQz 

M 


—x—M—l 


(3.25) 


with ^ = re ld , r > 0. For large M the integral has the leading form 


/ 71 

-1 


de z 


—n+M+x 


exp 


V«(*-))} = (-!) 


— n-\-M-\-x 


27rJ_ 


n-\-M+x 


(2 y/Q) 


uniform in n £ Z, where Jfc(x) denotes the Bessel function, and the result follows. 


□ 


Asymptotic formulas for {Rj(x)}, { < k ? (x')} appearing in the matrix elements of (12.1 hi) can 
now be obtained. 
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Corollary 2. Let M be even. We have 


Rm-2u(x + M) 


q=Q/M 2 


M/2—n 

(M -2n)\q-^ n+x ^ 2 ( ^ J 2 n + 2l +x ^VQ) 


1=0 


M/2-n 

-\f& ^2 •hn+2l-l+x{^\fQ) S j ■, 


(3.26) 


1=1 


R M - 2 n+i(x + M) ~ (M - 2n + l)!?-^*" 1 )/ 2 (j 2n _ 1+;c (2v / Q) - ^ 

q=Q/M 2 \ 


(3.27) 


$ M - 2 n(x + M) ~ a n (M — 2n)\q M ~ 2n ^ 2 V a v ' 2 J. U+X (2^Q), (3.28) 

q=Q/M 2 2 —' 


i/=—2n+l 


$M-2n+l(® + Af) 


q=Q/N 2 


a n+1/2 (M -2n + l)!<? M_2n+1)/2 

OO 

x ]T a "/ 2 (-J_„ +I (2\/Q) + 

Jv—2+x (2 a/Q))- (3.29) 

v=—2n 


Proof. The expansions CM and (1X771) follow immediately upon use of (1X71 in (1717) and 
(17X1 . In relation to {<!>_,■ (x)} we read off explicit formulas for the /3 n j as defined in (12.481) from 
(1717 . (1717 and substitute in (1717) to obtain 


$2n(®) 


< f > 2n+l (3?) 


r/2 V^L _ C 2 n)\q n (! ~ q) v+ W 

(1 - A /og) 2 7 2n + 1 (l - q) 2n ^ v\q v t 2 

v v=2n+l 


„x/2 _ 


a (2n + l)!g n +( 1 / 2 ) 


(1 - y/aq) 2 7 2n(l_ g )2n+l 


E 

v=2n 


(1 - q) v + 1 ~ 1 ' 


7 


u\q u / 2 


—C v {x) + 


(! ~<lf 


{v + 2 )(u + 1 )q 


C u+2 (x) 


(3.30) 


Use of (13.221) now gives (1X71 . (1X71 . 


Proposition 7. We have 



lim /9 (fe ) (hi + M,..., h k + M) 

M—> 00 


q=Q/M 2 


qdet 


S(hj,hi) I_(hj,hi ) 

D(hj,hi ) S(h h hj) \ jl=1 k 


□ 


(3.31) 


where 


S(x,y) = ^-(j x (2 v / Q)J y +i(2v / Q)-^(2v / Q)«/x+i(2 V / Q)) 

OO OO 

— u j/2 J-j+x(2\[Q) (</2/+2+j/(2-\/Q) — V«^ 2 /+i+ 2 /( 2 \/Q)^ ,(3.32) 

J=0 z=0 
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I(x,y) 


££« l/2 

n= 1 k =1 


Jn-\-x(^''J~Q') Jn— fc+j/(2 y/Q) Jn-\-y \/Q) Jn—k-\-x (2 y/Q))+e(x,y), (3.33) 


D(x, y) 

^ OO l 

= 7^EE( j ^(2v / Q) V^Wi^VQ)) ^2j+y-i (2t/Q) - ^/aJ2j+y{2^Q)^ 

v a i =i j=i 

-(x <-> y). (3.34) 


Proof. Using the results of Corollary 0 from the form of the matrix elements GI3-G2D 
we compute 


S(x,y) 


lim S(x + M, y + M) 

A-/—>oo 


q=Q/M 2 


and similarly for I(x,y), D(x,y). 


□ 


We note that Bessel function identities can be used to verify that the first term in (EUE1 can 
be written in a denominator free from 0 Prop. 2.9]. This allows (13.3211 to be replaced by 


S(x,y) = J2 J n +x (2^Q) J n+y (2^Q) 

n =1 

OO OO 

-^ai/ 2 J_ j+x (2^Q)Y^ (j 2 i+2+y(2v / Q) - V^zW 2 y/Q)) .(3.35) 
3=0 1=0 


□ 


The correlations dSTSU relate to the measure on diagrams of partitions (EH- a closely 
related measure, introduced in m , is 

P«>fl(A) = e -Vm-Q/2^Q^l /A) (3.36) 

where N = |A| and m = X^=i( —w ith denoting the length of the jth column of the 
diagram of A. This relates to decreasing subsequences in the involution. In the case (3 = 0 only 
rows of even length are permitted. For this model the corresponding fc-point correlations have 
recently been computed by Ferrari Lemma 5.2], in the context of the polynuclear growth 
model from a flat substrate. The results obtained have a very similar structure to that exhibited 
in Proposition 0 For example, the density (one-point correlation) is computed for (13.361) with 
(3 = 0 as 

/ 5 (?)’ 0) ^ = E( J «+-( 2 V / Q)) - J x+ i(2y/Q)(j2 4+2m-i(2v / Q) - 1 ~ 2 )> (3-37) 

n =1 m =1 

while for (13.711 with a = 0 Proposition 0 and (13.351) give 

oo 2 00 

= £ (^+*(2^)) -Jx&Vq) E Jx+2m (2 y/Q). (3.38) 

n= 1 m =1 

We also make mention that in the works Emnnsi the parameter dependent correlations have 
recently been calculated for other models generalizing the geometric picture of random involu¬ 
tions. 
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4 Asymptotic Correlation Functions 

4.1 De-Poissonization 

Following the pioneering work the following de-Poissonization lemma was derived in 0. 

Lemma 8. Let d > 0 be an arbitrary positive real number, and set 

= n + (2\/ d + 1 + 1) i/n log n, 
i/$f) = n — (2\/d + 1 + 1) log n. 

Let q = {g n i,n 2 }ni,ra 2 >0 be such that 


and define 


Qni,ri2 — 9711 + 1 , 712 ? 9711,712 — 9711,712 + 1 ? 0 — 9711,712 — ^ 

_ \711 \ 712 

_ —X 1 —X 2 „ A 1 A 2 


</>(Ai, A 2 ) = e 


dm 


n 2 I I 

ni!n2' 


ni,n2>0 

T/ierr f/iere exists constants C and uq such that for all ni,n 2 > no 

- cop + %'*) < «»,« < 4?) + cop + ( 4 . 1 ) 

In using this lemma, we take gn.m as the joint distribution of {Anln}j=i,...,l for random 
involutions of {1,..., N}, N = 2 n + m with n 2-cycles and m fixed points, and thus 

Qn,m = P r (A^m 5; a lj • • • t ^njm — a l)- 

We know from the argument of 0, Lemma 7.5] that this satisfies the monotonicity conditions 
required for the validity of Lemma [HI We further choose 

aj = 2 \fQ + Q 1 ^Sj, m = [Q/ 2], n 2 = [y/Q - 2wQ 1 ^ 3 ]. 

The lemma then gives 


lirn e~ Zl ~ Z2 

Q —>00 


~!L y'Ut 

J2 -T^T^^nin < < a l) 

z — J nlml 


n,m>0 


aj —2 Sj 

zi=Q/2 

Z2=y / Q—2wQ 1 / 3 


= lim Pr(A^) n < a u ..., < a t ) 

N —>00 


aj=2y/Q+Q 1 / 6 Sj ’ 

n=[Q/ 2] 

m =[\/Q — 2wQ 1 t 3 ] 


(4.2) 


assuming the limit on the left hand side exists. Regarding the left hand side, we know from 
El that it is equal to 


lim Pr (Q ’ Q) (Ai < ai,..., A; < a t ) 

Q —>00 


a. :j = 2 yfQ+Q 1 / & Sj 1 
yH=l-2ro/Q 1 / 6 


where Pr^ ,a ) refers to the ensemble defined by the 


measure 


E3). 


(4.3) 
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As demonstrated in (EH (EH, the probability in (EH is fully determined by the right hand 
side of font . Using font and fOTt this in turn can be written 

(a,Q) 


II i 1 -^ItrXXi-leir 


1=1 


r =1 
oo 


i + E 


(-i)” 


pi 


e n(E p { $) a) • • • > 


(4.4) 


p= 1 hi,...,/i p =0 i=l r=l 

Therefore our task is to compute the Q —■ ► oo limit of this quantity with 

I r = (a r ,a r - 1 ), a 0 = oo, aj = 2y/Q + Q 1/6 Si (j = 1,..., k), yfa = 1 - 2w/Q 1 ^ 6 . (4.5) 


4.2 The limit Q —> oo 

We will make use of a lemma due to Soshnikov IBS Lemma 2], in a form close to that used in 
0 . 

Proposition 9. Consider a sequence of point processes labelled by a parameter L. Suppose that 
for L —» oo, and after the linear scaling Xj i— > Ai,(xj) = + ai,xj of each of the coordinates, 

the sequence approaches a limit point process with correlations {p k }k=i, 2 ,... such that 



dx\ 


[ ^ n (e 

J — OO . 1 V _i 


i =1 r=1 


Xl, 


,Xk) 


o(k\), 


(4.6) 


and suppose furthermore that for each r = 1,..., l 

lim a k L dxi--- dx k p^AA^xf ),..., A L {x k )) = / dx i - / dx k p, k) (x x ,...,x k ). 

L ^°° JA L (Ir) JA L (I t ) V ’ Jlr Jlr 

(4.7) 

One then has 

i r =1 

OO / /*oo /*oo ^ ^ 

= 1 + E — / dx i - dx k Y[(^2^rX^)p(k)(xi,...,x k ). (4.8) 

fc=l • y -°° • y -°° i=l r=l 

The essential idea behind this result is that the condition (SB implies the moment problem 
for the number of particles in {I r } is definite, and thus the convergence of moments, which are 
integrals of correlation functions, implies the convergence of distributions. 

In using Proposition El the sequence of point processes will be those specified by the correla¬ 
tions (TUTTI) and thus labelled by the parameter Q. The limiting point process is that specified 
by the correlations appearing on the r.h.s. of EB- These correlations we will calculate to be 
(02). To deduce EH we must first verify EB for the correlations 031) . 

Lemma 10. There exists a positive number M such that 

p s $ ed (x i, ...,x k )< e —(*i+-+**) fcfc/2 M k (4.9) 

for all Xi E [so, oo) (i = were so is arbitrary but fixed. 
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Proof. For X —> oo we know that Ai(X) = 0(e 2x3/2 / 3 y Recalling (11.31) . we see from the 
explicit formulas EH that 


f 22 (X, Y) = f u (Y,X) = 0(e~ 2Y3l2/3 )0(e uX / 2 ), f 21 (X,Y) = 0(e uX/2 )0(e uY / 2 ), 

f 12 {X , Y) = o( e - 2x3/2 / 3 )o( e -2Y 3/2 /3). (4.10) 


The qdet can be expanded [HI 


& ( 0 ) 

qdet [/(Xj, Aj)]j J=1 = £ (_i)^n(/(A a ,X 6 )/(X 6 ,A c )---/(X d ,X a )) , (4.11) 

disjoint cycles 1 

eSk 

where the superscript (0) denotes the operation (l/2)Tr, and l denotes the number of disjoint 
cycles. We see from this that each term consists of factors of the form f u (X, X), f 22 (X,X), 
/ n (X,y)/ 22 (Y,X'), or f 12 (X,Y)f 21 (Y,X'). The bounds (TT~TT71) tell us that each term is 
bounded by 0(e~(' TlH h ' Tfc ^) (of course a sharper bound can be given, but this is sufficient 
for our purpose). 

We must also bound the dependence on k. For this we note that in general the replacement 
of the matrix elements of (11-131) by 

litv . . °P0 flit v v\ f!2 


r(x,Y) 

f 21 (X, Y) 


a{Y) 


f L1 (X,Y), / i2 (X,y) ~ a(X)a(Y)f l2 (X, Y) 


1 


a(X)a(Y) 


f 2L (X, Y) 


(4.12) 


leaves qdet unchanged. Choosing a(X) = e xl+ M (0 < fi <C 1) we see from (14.101) that each term 
is bounded for X,Y —► oo. According to Hadamard’s bound for determinants, the dependence 
on k is therefore bounded by k k / 2 M k for some M > 0. □ 

Because each I r is bounded from below, the bound (EH establishes (TPfll . The remaining 
task is to verify EZI)- First we make note of an alternative form of f 22 (X,Y). 


Lemma 11. The formula for f 22 (X,Y) in can be written 

DO 

Ai(s) ds 

roc v 

J e~ ut/2 Ai(t) dtj. 


f 22 (X,Y ) = K soft (X, y) + ^ ( A i(y) +1 


xe 


,uX/ 2(-u 3 /24 


(4.13) 


Proof. This follows by making use of the integral formula in EH for K soft , integration by 
parts, and use of the Fourier transform 


/ OO 

e~ vx Ai (x)dx = e~ y3/3 . 

-OO 


(4.14) 


□ 


It is precisely the form fljZEH which appears in the asymptotic form of S(x,y) relevant to 
verifying EH- 


Proposition 12. We have 
Q l / % S{2^Q + Q 1/e X, 2^Q + Q^Y) 


y/a=l—2w/Q 1 / G 


= r(x,Y) 


u=—Aw 


+ 0(Q~ 1 / 6 )0(e 


-Y\ 


(4.15) 
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Proof. Our main tool is the v —► oo asymptotic expansion 

J v (v-x{v/2) 1 ^) ~ f-) 1/3 Ai(x) + 0(-)0(e- a: ), (4.16) 

Vz// i/ 

uniform in x > xo (xo fixed). As noted in |6J Eq. (4.11)], this is a consequence of results due to 
Olver .5. • Using this with 

v = 2 v / Q + Q 1 / 6 X + n, x = X + nQ~ 1/e 


shows 


El </ti+x(2'\/Q)'4h-j/(2-\/Q) 

n=l 


x=2 v ^+Q 1 / 6 X 

y=2VQ+Q 1 / 6 Y 


-. OO 

= Ql73 E ( Ai ( X + n ^' 1/6 ) + 0(Q“ 1/6 )0(e- (X+nQ " 1/6) )) 

x (Ai(U + nQ~ 1/6 ) + 0(Q- 1/6 )0(e- (y+riQ " 1/6) )) 

= " Ai ( X + ^ Ai ( y + *) dt + 0(Q- 1/6 )0(e-(- Y+y ))), (4.17) 

where to obtain the final equality use has been made of the fact that a Riemann sum appears 
in the previous equality. 

The sum over j in the second term of (EP51) is not suited to the use of (Em To overcome 
this, we apply the generating function expansion 


oo , 

E t n J n {z) = exp(|(t--)) 


to conclude 


Y aj/2j -j+x^\[Q) = a x/2 _ 

l=o 

We can now use (n~T(fl) to deduce the asymptotic formula 


Y^ {j+x)/2 J j+x ^VQ))- 

3 = 1 


Ya j/2 J- j+x (zVQ) 


3= 0 


x=2y/Q+Q 1 / 6 X 
\/ct=l — 2w/Q^/6 


POO 

= e 8w3 / 3 e- 2wX 0 ( e - 2Xw2 /Q 1/6 ) _ / e 2wt Ai(X + t)dt + 0{Q- 1/6 )0(e~ x ). (4.18) 

Jo 

It remains to consider the sum over l in the final line of (13.3511 . For this we require in addition 
to (14.161) the uniform asymptotic expansion fjJ Below (4.12)] (again a consequence of results in 

123 ) 13 2 3 

J v±1 (y-x{yl 2) 1/3 ) ~ (^) 1/3 Ai(x)± Q 2/3 Ai \x) + 0^)0{e~ x ). (4.19) 

We find 


E ( J2l+2+y{^\fQ ) - y/aJ 2 l+l+y(2y/Q)^ 
1=0 


x=2VQ+Q 1 / 6 X 

y=2y/Q+Q 1 / 6 Y 

da=l—2w/Q 1 / e 


= 2^76 ( 2w + t ) dt ~ Ai (^) + 0(Q- 1/3 )0(e- y )). 


(4.20) 
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Substituting (HT71) . (ITTH1) and (14.2011 in (14.451) gives the stated result. □ 

For the remaining matrix elements in srm the following forms appear in the asymptotic 
form of D and I. 

Lemma 13. The formulas for f 12 and f 21 in irm can be written 

/ 12 (x,y) 


f\X,Y) 


JX JY 

(4.21) 

Proof. These follow straightforwardly upon using the integral formula in and the Fourier 
transform enD. □ 


1 /u d \ ru d 
4 V2 + ~dX) V2 + W 


x 


ds 


dt 


(Ai(F + s)Ai(X + t)~ Ai(X + s)Ai(y + t )), 


= -e“ |A '- y|/2 sgn(A - Y) 


a oo /*oo 

e u(Y-t)/2 A -soft( Xj dt _ J e u(X-t)/2 K soh^ t ) dt) 

roc /*oo 

_ e u(Y-X)/2 e -u 3 /24 / e us/2 Ai ^ ds + e u(X-Y)/2 e -u 3 /24 / 


Proposition 14. Let 0 < /x < 1/2 6e /ixed. We have 


q 1/3 d(2 V / q+q 1/6 x, 2 v^+g 176 ^) 


y/a=l—2w/Q 1 / e 


f U (X,Y) 


u=—Aw 


+0 


Q 1 / 6 


_ Vi+M —yi+M \ 

e ), 


7(2^ + q 1 /6 x, 2^q + g 1/6 y) 


= f 21 (X,Y) 


u=—Aw 



Proof. The procedure used in the proof of Proposition El suffices, with the error terms in 
gist ) and Aung) sharpened from 0(e A ") to 0 (e xl+>1 ). □ 

We remarked below the definition (EH) of a quaternion determinant that only specific 
combinations of the elements in the underlying 2x2 matrix can occur. We can use this fact 
together with the asymptotic expansions from Propositions 1121 and El to deduce an asymptotic 
formula relating p to p^ led - This in turn implies the Q —> 00 limit of (14.41) . 

Corollary 3. For Q —> 00 we have 


y/a=l—2w/Q 1 / 6 


Q k/6 p$ a) (2^Q + Q 1/e Xi, ... ,2 y/Q + Q^X k ) 

= P s $ ed (X u .. .,4 -Aw) + 0(^)0(e-( Al +-+^)). 

Consequently, with I r specified as in \4--5\) and so := 00 , 

(«,<?) 


(4.22) 


q 1 ™ (nti-y&xA, 


-leir 


1=1 


+E 


(-IF 


r =1 
r 00 


/a=l-2w/Q 1 / 6 
p k 


= 1 


P= 1 


dX i... r «,n(E . ,.,X r ; -4m). 

1=1 r=l 


(4.23) 
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Proof. Consider first Em>- The leading term is immediate. Regarding the error bound, as 
implied by the remark below Em, the elements in the 2x2 matrix of appear as factors 

in on only in the combinations S(x, x ), S(x, y)S(y, x') or I(x, y)D(y , x'). This, together with 
the individual error bounds in Propositions 1121 and 1141 imply the error bound in (14.221) . 

We know from Proposition El that to establish (14.241) it suffices to establish (14. (il) and Oil . 
As previously remarked, the bound (EH implies EH. Thus it remains to check that with cij 
as in (14.51) . 


Jin^ Q fc / 6 ]T H X^a^/S) fa ’''' ’ ' h < 

h\,...,hp=0 i=l 

V 


yja=l-2w/Q V 6 


/ °° /*00 r 

dx 1 • • • / dx p JJ • • • > X V 1 ~ 4w )- ( 4 - 24 ) 

-OO J —O O i=\ 


For this we note 

OO p 


S II X (2y/^+Q 1 ^s r ,2y/^+Q 1 ^s r -i) P (p) fa ’ ' ' ■ ’ fa 

hi,...,h p =0 i= 1 

(i) (Q,a), 


£ ( 2 v^+« 1 '‘ i -- 2 v , e+« 1/t = 

*1> —•*?: i= 1 

2v'^+Q 1 /%-SZ 

Recognising this as a Riemann sum, and substituting EiH gives Em. 


□. 


4.3 Proof of theorem 1 

According to EH , EH , (EH and EH 


lim Pr(Ag^ < a\,..., X n,m < 


TV —>00 


a,j — 

n=[Q/2] 

ra= [v^ - 2it)Q 1 / 3 ] 


= lim Pi 4Q ’ a) (Ai < ai,..., A* < a t ) 
Q —>OG 


dj=2y/Q+Q^-/^Sj 
y/ac=l — 2w/Q^/® 


y 


(l)Er=l n r QT, l j=l n j 




nil ■■■nil d£ ni ■ ■ ■ d£ ni Q^oo . 

(ni ,...,n;)eL; J-l r —1 


(«,Q) 


a=l—2UI/Q 1 / 6 

(4.25) 


The limit in the final expression is evaluated according to Em , which shows Em can be 
written 

]T E^dnr^SrySr^jr^.P, (4.26) 

where E( w \{n r , (s r , s r - i)}r=i,...,z) denotes the probability that the interval (s r ,s r _i) contains 
no particles in the point process specified by the correlations (11.141) with w given by (11.111) . But 
Em is, according to EH , equal to the distribution function F 0 (si,..., si; w) appearing on 
the right hand side of 031) . 
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